Influence of magnetic viscosity on domain wall dynamics under spin-polarized currents 
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We present a theoretical study of the influence of magnetic viscosity on current-driven domain 
wall dynamics. In particular we examine how domain wall depinning transitions, driven by thermal 
activation, are influenced by the adiabatic and nonadiabatic spin-torques. We find the Arrhenius 
law that describes the transition rate for activation over a single energy barrier remains applicable 
under currents but with a current-dependent barrier height. We show that the effective energy 
barrier is dominated by a linear current dependence under usual experimental conditions, with a 
variation that depends only on the nonadiabatic spin torque coefficient j3. 

PACS numbers: 75.60.Ch, 75.75.+a, 85.75.-d, 05.40.-a 
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I. INTRODUCTION 

Spin-polarized currents flowing through ferromagnetic 
media exert torques on local magnetic moments in re- 
gions where a spatial magnetization gradient occurs.— 
This interaction originates from an exchange interaction 
between the charge carrier spins and local magnetiza- 
tion, and has been observed in systems based on tran- 
sition ferromagnet metals 2 -^! 4 - and dilute magnetic semi- 
conductors^ alike. In the context of domain walls, such 
spin-torques can be sufficiently large to induce wall mo- 
tion even in the absence of magnetic fields £ Interest in 
current-driven domain wall dynamics stems from both 
fundamental considerations, whereby the problem in- 
volves reconciling complex spin-dependent transport pro- 
cesses with nonlinear magnetization dynamics, and po- 
tential applications for emerging spintronic technologies, 
such as novel logic circuits 7 - and nonvolatile memory de- 
vices jSis 

The salient features of current-driven domain wall dy- 
namics are captured in the modified Landau-Lifshitz 
equation of motion for the unit magnetization m ; 10 ' 11 

dm -j _ dm _ ^ 

— +7omx H c g — amx -—- — — (u- V)m+pmx (u- V)m, 

where the right-hand side describes the current- 
dependent terms. Here, u represents an effective 
spin current drift velocity (see Section II for defini- 
tion). The first term on the right hand side de- 
scribes the adiabatic contribution, which arises from 
transport processes in which the charge carrier spins 
track the local magnetization variations as they traverse 
the wall. In contrast the second term, proportional 
to (3, represents nonadiabatic effects. There remains 
much debate over the magnitude of the nonadiabatic 
torques, 12 ! 13 ' 14 ! 15 ! 16 ! 17 ' 18 ! 19 ! 20 ! 21 ! 22 ' 23 whereby theoreti- 
cal estimates invoking momentum transfer ) 12 ' 21 i 22 spin- 
mistracking , 13 ' 17 spin-flip scattering , 16 ' 18 ! 20 ' 22 or even ap- 
parent nonadiabatic effects involving magnon emission, 23 
differ at least by an order of magnitude. The nonadia- 
batic term is important because it allows a domain wall to 
be driven into the viscous regime^solely by applied cur- 



rents with no current thresholds in perfect systemsj 11 ' 24 
This is in stark contrast to the effect of a pure adia- 
batic torque, whereby a current threshold governed by 
the hard-axis anisotropy needs to be surmounted before 
pure current-driven wall motion is possible.— Such dy- 
namics can be readily deduced from the reduced equa- 
tions of motion for a one-dimensional Bloch wall, which 
can be obtained from (1) by assuming a specified wall 
profile and integrating over the spatial coordinates; 1 - 1 - 

Xq/\ — a<j> — r ynoH± sin<^>cos</> + u/X, (2a) 
+ aX o /X = -ffi o H a +0u/X, (2b) 

where Xo denotes the wall position, <p describes a tilt an- 
gle that is associated with the conjugate momentum to 
Xq, X is the domain wall width, and H±_ is a transverse 
anisotropy field. From these equations it is immediately 
apparent that the nonadiabatic term plays a similar role 
to that of an applied magnetic field H a , which has led 
to some attempts of quantifying (3 experimentally by as- 
sociating field-like variations, such as a reduction in the 
propagation field under currents, directly with the nona- 
diabatic torque. 

Despite the lack of consensus over the theoretical mag- 
nitude of the nonadiabatic torque, it is generally expected 
that large nonadiabatic effects should appear in systems 
with narrow domain walls. Good candidates to exhibit 
large nonadiabatic effects are therefore ferromagnetic al- 
loy films with perpendicular anisotropy, such as Co/Pt 
and Co/Ni multilayer a 25 ! 26 ' 27 or FePt alloys^ 8 - as wall 
widths in such materials can reach down to near-atomic 
dimensions of 1 run. Furthermore, wall structures in 
these materials are expected to closely resemble Bloch 
domain walls, because the shape anisotropy associated 
with magnetization perpendicular to the film plane has 
the same symmetry as the intrinsic uniaxial anisotropy. 

However, it is well established that wall motion in per- 
pendicular anisotropy materials is largely dominated by 
intrinsic defects, to the extent that pinning fields can po- 
tentially be larger than the Walker breakdown fieldi^ 9 - For 
this reason, it is expected that thermal activatio n 30 ! 31 
and magnetic viscosity^ 2 , to be important for current- 
driven dynamics in such systems. In this context, the 
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basis for the association of the nonadiabatic term with 
an effective field is unclear. While for the streaming dy- 
namics the P term enters the equations of motion like a 
field, a linear relationship between an applied magnetic 
field and an energy barrier, associated with pinning po- 
tentials due to defects, for example, does not necessarily 
exist in general.— In this article we show that an Arrhc- 
nius law remains valid for describing transition rates r _1 
associated with thermally activated wall processes under 
spin-polarized currents, r = To exp^E^/hsT), but with a 
current-dependent effective energy barrier that depends 
only on nonadiabatic torques, Ef, — Eb(j3I). Under usual 
experimental conditions, we show that the dominant con- 
tribution is a linear current term, which suggests that a 
direct experimental measure of f3 can be obtained from a 
linear variation of the characteristic Arrhenius rate with 
current for domain wall depinning from a defect. As a 
corollary, the association of (3 with an effective magnetic 
field is only possible if the energy barrier varies linearly 
with applied magnetic fields. 

This article is organized as follows. In Section II, 
we present the one-dimensional Bloch wall model for 
wall dynamics. In Section III, we derive the Langevin 
and Fokker-Planck equations for the stochastic wall dy- 
namics, and solve for the stationary probability dis- 
tribution function. In Section IV, we calculate the 
current-dependent Arrhenius rate equation for thermal- 
activation of the domain wall over an energy barrier, and 
apply these results to depinning processes associated with 
point defects in Section V. We present some discussion 
and concluding remarks in Section VI. 



II. ONE-DIMENSIONAL BLOCH WALL MODEL 

We consider the effects of magnetic viscosity on the 
dynamics of a one-dimensional Bloch wall under spin- 
polarized currents. Bloch wall structures are expected 
to arise in magnetic films with perpendicular anisotropy, 
as the stray dipolar fields in this geometry act only to 
renormalize the unaxial perpendicular anisotropy. We 
will concern ourselves with the one-dimensional motion 
of a Bloch wall in a magnetic wire. Let z represent the 
axis perpendicular to the film plane, which is collinear 
with the axis of uniaxial anisotropy K u . Let X rep- 
resent the direction of the long axis of the wire, which 
corresponds to the direction in which the magnetization 
variation constituting the domain wall arises. It is as- 
sumed that the magnetization remains uniform along the 
y and z directions along the wire. Using spherical coordi- 
nates to represent the spatial magnetization distribution 
m(X) = (mx,'my,m z ) = (cos </>sin#, sin0sin 9, cos#), 
where [9, <fi] — [9(X), <f>(X)], the simplest magnetic 
Hamiltonian leading to a domain wall is 



EL, = dV A 



drrii drrii 
~dX~dX 



K u m: 



(3) 



where the first term represents the exchange interaction 
(with summation over repeated indices implied) with the 
exchange constant A, and the second term represents the 
uniaxial anisotropy. The static Bloch wall profile can 
be obtained using a variational procedure by minimiz- 
ing the energy functional (J3|) with respect to [#,</>]; one 
readily finds 9(X) — 2 tan -1 exp[— (X — Xq)/X], where 
the characteristic wall width is A = \J Aj K u . Note 
that in magnetic films with perpendicular anisotropy, 
the constant K u represents an effective anisotropy K u — 
K u .o — that represents the difference between 

the magnetocrystalline term K u ^ and the demagnetiz- 
ing fields due to the film geometry. 

The salient features of the magnetization dynamics of a 
domain wall can be accounted for by the one-dimensional 
Bloch wall model. In this approach, the dynamics is sim- 
plified by integrating out the spatial degrees of freedom in 
the Landau-Lifshitz-Gilbert equation of motion (fl]) . This 
is achieved by elevating the domain wall position Xo(t) 
to a dynamical variable and by assuming that (f> — <p(t) is 
spatially uniform. Within this approximation, the wall 
dynamics is entirely parametrized in terms of the two 
wall coordinates (Xo,4>). In the absence of thermal fluc- 
tuations, the equations of motion are given by (including 
the adiabatic and nonadiabatic terms),— 



I dX 90 

x dt a dt 

def) a dX 
~dt + ~X~dT 
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where u = jPg[iB/2eAd s represents an effective spin cur- 
rent drift velocity, which accounts for the spin-transfer 
torque, 7 is the gyromagnetic constant, and S c is the 
cross-sectional area of the wire. P is the spin polariza- 
tion, g is the gyromagnetic ratio, /is is the Bohr mag- 
neton, e is the electron charge and M s is the saturation 
magnetization. The energy E(Xq, <f) represents the total 
magnetic energy of the domain wall, which includes the 
Zeeman energy from applied and stray dipole fields, and 
contributions from local field and anisotropy variations 
due to defects. 

This set of coupled differential equations can be simpli- 
fied using dimensional analysis by scaling out the charac- 
teristic length, time, and energy scales. The domain wall 
width A is the obvious characteristic length scale, which 
leads us to define the reduced spatial variable x 



Xo 
A ' 



(5) 



A dimcnsionlcss time variable can be defined in terms of 
the characteristic frequency w m — 7/ioM s , 



t — U! m t. 



(6) 



In a similar way, we can define a dimensionless energy 
e = E/Eq, where the choice of 
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allows the equations of motion to be simplified neatly. 
For a ferromagnetic material with HqM s — 0.75 T and 
A = 10 nm, patterned into a wire geometry with a cross- 
sectional area of S c = 100 nm x 1 nm, Eq corresponds to 
a temperature of 6.5 x 10 4 K. In the same spirit, we will 
find it convenient to define a dimensionless temperature 
f = T/Tq, where T = E /k B . 

With these definitions and u = u/(\u> rn ) being the 
dimensionless spin drift velocity, we obtain the dimen- 
sionless form of the equations of motion, 



dx 

dv 
dt 



dv 
dx 

l Tt 



de 
dv 

OX 



(8a) 
(8b) 



where we have used the fact that the wall angle </> is 
linear in wall velocity at low velocities, which allows us 
to associate <j) with a dimensionless velocity v, v = <f>. 
This association makes the analogy between the ID Bloch 
wall dynamics and the kinematics of a point particle more 
explicit. In this context, we note that the kinetic energy 
of the domain wall is associated with the excursion angle 
4> out of the wall plane thorugh the hard-axis anisotropy, 

Eun = ^aM'j J dV sin 2 8 sin 2 = fj, M*S c \ sin 2 0, 

(9) 

which upon substituting the characteristic energy scale 
and u, and including a particle potential U(x), leads to 
the reduced domain wall energy 



e = \v 2 + U(x) 



(10) 



which is valid at low velocities (sin 2 cf> w </> 2 = v 2 ). The 
term U (x) contains the Zeeman energy and any pinning 
potentials. 



III. LANGEVIN EQUATIONS AND 
FOKKER-PLANCK THEORY 

In order to describe processes involving thermal acti- 
vation, we need to extend the theory to include thermal 
fluctuations. This can be achieved by adding stochastic 
forces rj to the equations of motion^ 



dx 

~dt 

dv 

dt 



a 



dv 

"at 

dx 



de 
dv 



de_ 

dx 



u + r) v (t), 



(11a) 
(lib) 



where the rj terms represent white Gaussian processes 
with zero mean, {rj{t)) = 0, with two-time correlation 
functions of the form (^(i)^^')) = (2aT)(J /XiV 5(t — 
<'). This definition of the noise is consistent with the 
fluctuation-dissipation theorem for the closed thermody- 
namic system u — 0. 

Next, we seek to derive the Fokker-Planck equation 
that is associated with the coupled Langevin equations 



in Eq.[TTJ The Fokker-Planck equation describes the time 
evolution of the probability density P(x,v,t) of finding 
the system in a state (x, v) at time t. We proceed by first 
rewriting (|11[) to read 



(1 + a 2 ) 
(l + o, 2 ) 



dv 
dt 

d.r 



a— ) + u((3 — a) 



de 

^dx dv 
+(rj x -ar) v ), (12a) 

dv dx J 
+(ar) x + r) v ). (12b) 



Now, for a set of coupled nonlinear Langevin equations 
of the form, 



±i = hi({xi},t) + g lj {{x i },t)r]j{t), 



(13) 



where Xi represent the state variables and rji are the ran- 
dom fields, the corresponding Fokker-Planck equation for 
the probability density P of finding the system in a state 
{xi} at time t, 



dP{{x,},t) 
dt 



= L FF P{{xi},t), 



(14) 



with the Fokker-Planck operator 
Lpp - 



d d" 2 
— D l ({x l },t) + ^-^-D t3 ({x z },t), (15) 



can be constructed from the drift coefficients 

d 

Di({xi},t) = hi({xi},t) + gk]{{xi},t)-^—g ij ({x i },t) 



and diffusion tensors 

Dij({xi},t) = g ik ({x i },t)g : j k ({xi},t), 



(16) 



(17) 



where summation over repeated indices is implied^ By 
applying this prescription to the Langevin equations in 
(I12p . we find for the domain wall system 



Lpp — Lq + L u , 



(18) 



where 
L 

Li, 



d f de de 
dx \ dv dx 



d_ (de_ 

dv \ dx 



de 
dv 



& 



2 



d 2 



dx 2 dv 2 



(19a) 
(19b) 



which agrees with the one-dimensional low-velocity limit 
of the operator found by Duine et ali^ Without loss of 
generality, we have assumed a 2 <C 1 to simplify the no- 
tation. 
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In the absence of currents, it can be easily verified that 
the Boltzmann distribution function satisfies the station- 
ary Fokker-Planck equation, 



L P (x,v) = 0, 



where 



P (x,v) = Z 



(20) 



(21) 



with Z being a normalization constant. In the presence 
of currents, it is also possible to find the stationary solu- 
tion to the full Fokker-Planck equation. From inspection 
of (jlip . we note that the adiabatic torque can be assim- 
ilated as a shift in the wall energy e — > e + uv, while 
the nonadiabatic torque leads to the shift e — * e — uj3x, 
leading to the ansatz 



Po(x,v) 



exp 



e + u(—/3x + v) 
T 



(22) 



which is readily verified to be a solution to the stationary 
Fokker-Planck problem for u ^ 0. 

This distribution function shows that the Gaussian dis- 
tribution over the particle "velocity" v becomes centered 
around the value v = —u as a result of the adiabatic 
torque, Pq(v) oc exp[i(t> + u) 2 ]. This would appear to 
be an inconsistent result, whereby the net domain wall 
flow takes place in a direction opposite to the spin cur- 
rent. The origin of this apparent conundrum lies in the 
association of sin(</>) with a wall velocity v. It is well 
established that while the adiabatic torque cannot drive 
the wall into a streaming motion at steady state (such as 
that achieved with a magnetic field or with a sufficiently 
large adiabatic torque^), it does lead to a transient dis- 
placement of the wall with a finite value of cj> at rest. 
The shift in the Gaussian distribution is a manifestation 
of this phenomenon, which restates that the tilt angle <fi 
attains a finite value under adiabatic torques, but does 
not correspond directly to a finite streaming wall velocity. 



IV. ARRHENIUS TRANSITION RATE 

To calculate the thermally-activated depinning transi- 
tion rate for a domain wall pinned at a defect, we follow 
the transition rate theory due to Kramers i 36 ' 37 This ap- 
proach relies on drawing a strong analogy between the 
stochastic dynamics of the domain wall with the Brown- 
ian motion of a point particle in one-dimension. We con- 
sider a particle in a potential well as shown in Figure [TJ 
with an energy barrier given by e^o = U(xb) — U(x a ) > 0. 
The transition rate r _1 of particles out of the potential 
well can be estimated by^ 7 - 



Jb_ 



(23) 



where jb is the particle current at the summit of the 
potential barrier x = Xb, 



Jb 



dv J x (x b ,v), 



(24) 




FIG. f: (Color online) Escape of a particle over an energy 
barrier by thermal activation. The barrier height is = 
U(xb) — U (x a ) > and the characteristic potential well width 
is Ax = Xb — x a > 0. 



with J x being the probability current density, and n a is 
the particle population in the well, 



;/,, = / dx 

'well 



dv Pq(x, v). 



(25) 



To compute the well population, it is convenient to 
use a quadratic expansion of the potential well around 



U(x) 



U(x a ) + -L0 2 a {x- 



Xaf 



(26) 



where to 2 > is the square of a characteristic angular fre- 
quency that describes the curvature of the potential well 
at x — x a - We assume that the particles injected into the 
well possess energies at least a few ksT below the barrier 
height and are rapidly thermalized before thermal acti- 
vation processes take place, which allows us to use the 
steady-state distribution to describe the particle popula- 
tion at the bottom of the well. By using the stationary 
distribution function given in (|22p . we find 



1 2ttT 

Z LO a 



(27) 



x exp <; j 



-U{x a ) +u/3x a + — 1 + -^ 



To describe the particle current at the barrier summit, 
we follow Kramers' approach by seeking a modified so- 
lution to the stationary Fokker-Planck equation of the 
form 



P (x,v) =Z- 1 C(x,v) 



exp 



e + u(—j3x + v) 



T 



(28) 



where the function £(x, v) accounts for sources, i.e. the 
injection of particles into the well at x = x a , and sinks, 
whereby particles leaving the well x > Xb are removed 
from the system. As before, we can approximate the 
potential energy at the barrier summit as 

U{x)^U{xb)~\u 2 b {x-Xb) 2 , (29) 

where u> 2 > 0. The limiting behavior required for £(&, v) 
is that it should be I in the well and tend toward zero for 



5 



x > xi,. This behavior is only possible if £(x, u) depends 
on a linear combination of x and v only. 3 ® Let 

£ = (x — x b ) + a(v + u) H J • 



where we have defined 



^6 = 



2u b 



a{l+u 2 ) 



(37) 



By substituting this solution ansatz into the stationary 
Fokker-Planck equation, we find a simple second order 
differential equation for Q in terms of £, 



d 2 t dC 



0, 



where 



A+ = 



(a - a±)uj 2 
(1 + a 2 ±)af 



(30) 



(31) 



and a± correspond to the roots of the quadratic equation 
1 + aa± — a±u> 2 (a± — a), 

The physically correct solution that satisfies the appro- 
priate boundary conditions for £ requires a~, giving 



C(0 = ^rfc (v^Ce) 



(33) 



where erfc(a;) is the complimentary error function, 

ev{c(x) = (2/^)Cdte- t2 . 

The probability current density describing the flow in 
the spatial coordinate x is defined in terms of the drift 
and diffusion coefficients as 35 



d 

J x — D X P — — D X jP. 

UXj 



(34) 



By inserting the stationary Fokker-Planck solution sub- 
ject to the boundary conditions defined by the function 
C(£), we find the overall particle current to be 



3b 



-aT 



dv (v + u)Po(xb, v 



dv^le-^. (35) 



By using the dummy integration variable z = v + u and 
the identities 



dz ze az erfc(6z + c 

dz e -z 2 e -Kz+cf = 



be ~ac 2 /(a+b 2 ) 



tVa + b 2 

j-abc 2 /(a+b) 



we find an analytic expression for the particle current at 

Xb, 



1 



Jb = 2z exp 1r 



i 



u 2 ( f3 2 

-U (x b ) + uf3x b + — 1 5 

2 V uj 2 



x-(«t [i-n b -Lu 2 b (i + n 2 )]), 

UJb 



(36) 



to simplify the notation. 

After combining (|28p and (|3"6"|) with ([2^)1 . we recover an 
Arrhenius law for thermal activation over a single energy 
barrier, 



f 1 

— = — exp 
r to 



f 



(38) 



but with a current-dependent effective energy barrier 



e b (u) = e 6! o - {f3Ax)u + 



1 ( 1 



f 



mf (39) 



with linear and quadratic terms in the applied current 
(through u) , and an attempt frequency of the form 



f 

to 



2tt 



[i-r> fc -^ b 2 (» b -i)] 

2uj b 



(40) 



Eqs. I3"81l4"01 constitute the main result of this paper. Be- 
cause u <C 1 for typical low-current measurements in 
experiment, the dominant contribution to the current- 
dependence of the energy barrier is expected to be the 
linear term. It is important to note that both the linear 
and quadratic terms depend only on the nonadiabatic 
coefficient f3, which suggests that nonadiabatic torques 
can be quantified experimentally by measuring how the 
energy barrier varies with applied currents. 



V. APPLICATION TO POINT DEFECTS 

In this section, we apply the theory developed in the 
preceding section to a simple model of domain wall pin- 
ning. Two kinds of point defects are considered: (1) A 
reduction in the uniaxial anisotropy. (2) A local hard axis 
pinning field. Both defects lead to a local minimum in the 
potential energy of the domain wall, which leads to pin- 
ning of the wall at the defect. From this simple model, it 
is possible to calculate directly important parameters for 
the Arrhenius transition rate, such as the field and cur- 
rent variation of the attempt frequency, barrier height, 
and barrier curvatures, which facilitates comparison of 
the present theory with experimental data. 

For the anisotropy defect, which is characterized by a 
pinning strength Hp d and assumed to be located at the 
origin, the corresponding energy is 



E. 



ad 



z adl x H; d M s / dV sm 2 0[{X-X o )/X]S(X), 



adt xoH^M s S c X sech z (x). 



(41) 



Despite the pointlike nature of the defect the potential 
well seen by the domain wall has a finite width because 
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FIG. 2: (Color online) (a) Domain-wall potential energy U 
for three different applied fields with H p /M s = I. (b) Field 
dependence of the energy barrier e^o, where squares corre- 
spond to the numerically calculated values and the solid line 
represents a fit of the form e^o = C(l — fe M )". 



the wall has a finite spatial extension^ For the field 
defect, which is characterized by H^, a similar expression 
is obtained, 



E 



I'd 



z U fx H; a M s / dV sm9[{X - X )/X}d(X), 



z id fJ, Hl d M s S c X sech(x). 



(42) 



In the presence of an applied field H a along the positive z 
axis, there is an additional Zeeman term in the potential 
energy of the form 



Eg — -2fi H a M s S c X . 



(43) 



The coefficient z in each of the pinning energies is chosen 
such that the energy barrier vanishes for H a = H p . In 
terms of the scaled field h = H a /H p , the total (scaled) 
potential energy U — (l/Eo)(Ed + E z ) for the domain 
wall are 



and 



U ld (x) 



^(^sech 2 (z)+te) (44) 



— 2- (2sech(x) + hx) . 
M. 



(45) 



for the anisotropy and field defects, respectively. 

The potential profile for both anisotropy and field de- 
fects is shown in Fig. [3^ andHt, respectively, as the ap- 
plied field is varied from zero to a magnitude equal to 
the pinning field at which the barrier disappears. It is 
straightforward to calculate the field dependence of the 
energy barrier with these pinning potentials. These re- 
sults are shown in Fig. [^b anddJi. The field variation of 
the energy barrier can be described empirically by 
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e 6 . = C(l - h»y 



(46) 



FIG. 3: (Color online) Field depedence of (a) barrier width 
Aa; and (b) characteristic frequencies u) a ,b- The latter are 
proportional to the curvature of the potential well and en- 
ergy barrier, respectively. Solid lines represent the anisotropy 
defect, dashed lines represent the field defect. 



where li, v are exponents. From fits using this equa- 
tion, we note that v ~ 3/2, which corresponds to a 
strong pinning form. 33 For each defect, we have also in- 
cluded the simple linear field variation for the energy bar- 
rier— for the purposes of comparison. While the linear 
form captures the general trends qualitatively, there is a 
large quantitative discrepancy between the exact form of 
£b,o(H) and the linear approximation, at least for the 
point defects considered here. Point defects therefore 
constitute an important class of pinning potentials for 
which a linear field/current equivalence does not apply 
for nonadiabatic torques. 

Using this model, it is also possible to easily determine 
the field variation of the width of the potential well and 
the curvature of the potential well and barrier summit, 
given by the frequencies u> a and u>b, respectively. These 
are presented in Fig. [3J The width of the barrier is of the 
same order of magnitude as the wall width A, as expected, 
with a weak dependence on h over a large range of ap- 
plied fields. Because the linear current term in the energy 
barrier is also proportional to Ax, it may be important 
to account for such field variations in comparisons with 
experimental data. Nevertheless, one can still obtain a 
good order-of-magnitude estimate of nonadiabatic effects 
by assuming Ax ~ I. 

Some of the key features of the influence of spin torques 
on the energy barrier (|39[) are presented in Fig. U for an 
anisotropy defect. In Fig. andSh, the field variation 
of the energy barrier is shown for positive and negative 
currents, respectively, for different values of the nona- 
diabatic torque coefficient. While the variations to the 
zero current limit are relatively small over most of the 
field range, there is a significant contribution from the 
quadratic current term near h = and h = 1 at which 
the curvature of the barrier summit u>b tends toward zero. 
This behavior is particularly explicit for the large /3 = 1 
case considered. We reiterate that the present theory is 
only valid regions in which a well-defined potential well 
exists, i.e. for uj a ^ ^ 0; the divergences at h = and 
h = I are therefore unphysical. Nevertheless, we do ex- 
pect that the quadratic current term to become impor- 
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FIG. 4: (Color online) Energy barrier tb (in reduced units) 
as a function of (a,b) applied field h, with (a) u = 0.1 and 
(b) u — —0.1, and (c,d) applied current u, with (c) h — 0.1 
and (d) h — 0.5, for three values of the nonadiabatic torque 
coefficient j3. 



tant, particularly for large j3 values, both at small fields 
and at fields close to the depinning field. This can be 
seen more explicitly in Figs. and d in which the cur- 
rent dependence of the energy barrier is shown at weak 
and intermediate fields, respectively. 



VI. DISCUSSION AND CONCLUDING 
REMARKS 

We have found that the transition rate associated with 
domain-wall depinning transitions driven by thermal ac- 
tivation, under applied magnetic fields and spin-polarized 
currents, can be described by an Arrhcnius law in which 
the effective energy barrier depends on applied currents 
/ in the following way, 

E„(I) =E b , -a 1 I + (T 2 I 2 , (47) 

where the two efficiency parameters are defined as 

ai = (3P-^f, (48a) 
e A 

Under usual experimental conditions the dominant con- 
tribution is expected to come from the linear term, which 
is proportional to the nonadiabatic coefficient (3 and the 
spin-polarization of the material. The linear term also 
depends on the ratio between the width of the potential 
well and the domain wall width, AX/A, but this ratio 
is expected to be of the order of unity for pointlike de- 
fects. The quadratic current contribution can become 
important for small applied fields and for fields close to 
the switching field, where the curvature at the energy 



barrier summit (u>%) becomes small. However, we ex- 
pect that it would be unrealistic to measure the statis- 
tics of depinning under vanishing magnetic fields (i.e., 
large energy barriers) as the transition rates under these 
conditions would be immeasurable in practice. In com- 
parison with the effective temperature model proposed 
for thermally-assisted magnetization reversal with spin- 
transfer torques in magnetic nanopillars^ we note that 
spin-torque effects here appear also to renormalize the 
energy barrier (or the effective temperature), but with 
a different functional form that cannot be expressed in 
terms of a critical current (1 — I/I c )- 

The transition rate has been obtained by assuming low 
domain-wall velocities for which we have assumed that 
the tilt angle <j> is proportional to a streaming wall veloc- 
ity v. This approximation is valid for wall propagation 
below Walker breakdown^ above the Walker threshold, 
wall motion is oscillatory by which the angle <f> under- 
goes a precessional motion. In metallic alloys with per- 
pendicular anisotropy, it is possible that pinning fields 
exceed the Walker field ^ which means that wall propa- 
gation would immediately proceed in the Walker regime 
following depinning from a defect. We believe that the 
present theory remains valid even for large pinning fields 
exceeding the Walker field, because the motion associated 
with the depinning transition is more likely to be viscous 
rather than precessional. This can be understood as fol- 
lows. The Walker regime is attained when the slope of the 
domain wall potential —dU/dx, which corresponds to a 
linear force on the wall, exceeds a certain value. Because 
the slope of the potential well decreases to zero at the 
pinning center, any domain wall motion driven by ther- 
mal fluctuations around this position is likely to result in 
viscous motion rather than precessional motion because 
the net force acting on the wall remains small. This is 
equally true at the barrier summit. Because the tran- 
sition rate is ultimately determined by the dynamics at 
the bottom of the potential well and at the barrier sum- 
mit, we expect the approximation of viscous motion will 
remain relatively robust in spite of possible thermally- 
driven excursions into the Walker regime. 

We wish to reiterate that the current-induced modifica- 
tions to the energy barrier do not appear necessarily as an 
effective field. This is not entirely surprising given that 
the nonadiabatic torques enter only as a fieldlike contri- 
bution in the equations of motion, but cannot themselves 
be derived from a potential because they are nonconser- 
vative torques. As such, we contend that measures of the 
(3 coefficient through fieldlike variations, e.g. associating 
a change in propagation field with a nonadiabatic torque, 
can lead to erroneous estimates if thermal activation or 
magnetic viscosity is not accounted for properly. This 
problem is particularly important for ferromagnetic al- 
loys with perpendicular anisotropy where we expect the 
influence of magnetic viscosity, due to the large contri- 
bution of intrinsic defects, to be important for current- 
driven wall dynamics. We suggest that measurements 
of the current-dependence transition rate of domain-wall 
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depinning from defects should provide a more direct and 
accurate determination of the nonadiabatic torque in ex- 
periment. 
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